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Abstract 

Mycorrhizal association is known to be innportant to orcllid species, and a complete understanding of the fungi tliat form 
mycorrhizas is required for orchid ecology and conservation. Liparis japonica (Orchidaceae) is a widespread terrestrial 
photosynthetic orchid in Northeast China. Previously, we found the genetic diversity of this species has been reduced recent 
years due to habitat destruction and fragmentation, but little was known about the relationship between this orchid species 
and the mycorrhizal fungi. The R/?/zocton/a-like fungi are the commonly accepted mycorrhizal fungi associated with orchids. 
In this study, the distribution, diversity and specificity of culturable /?h/zocton/o-like fungi associated with L. japonica species 
were investigated from seven populations in Northeast China. Among the 201 endophytic fungal isolates obtained, 86 
/?/?/zocfor)/o-like fungi were identified based on morphological characters and molecular methods, and the ITS sequences 
and phylogenetic analysis revealed that all these Rh/zocfon/o-like fungi fell in the same main clade and were closely related 
to those of Tulasnella calospora species group. These findings indicated the high mycorrhizal specificity existed in L japonica 
species regardless of habitats at least in Northeast China. Our results also supported the wide distribution of this fungal 
partner, and implied that the decline of L. japonica in Northeast China did not result from high mycorrhizal specificity. Using 
culture-dependent technology, these mycorrhizal fungal isolates might be important sources for the further utilizing in 
orchids conservation. 
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introduction 

The mycorrhizal association is ubiquitous but very important 
symbiosis in nature, which plays an essential role in the 
maintenance of most terrestrial ecosystems [1]. Over 90% of all 
plant species can form mycorrhizas with different kinds of fungi, 
and the existence of mycorrhizal fungi can confer to their hosts 
many adaptive advantages via improved water and nutrient/ 
minerals uptake from the soil [2,3,4], enhanced plant growth [5,6], 
reduced toxic element accumulation [1,7], and increased resis- 
tence to pathogen damage [8] . Nowadays, based on their essential 
role in plant species and the whole ecosystems, more attention has 
been given to the potential of exploiting fungi for effective use 
[9, 1 0] . Although a large number of mycorrhizal fungi existed in 
plants, only a fraction has been described and explored to date, 
and their function in ecological systems remains indistinct. 

The Orchidaceae, which is one of the largest and most diverse 
plant families, is distributed worldwide [11]. However, many 
orchid species have suffered dramatic declines in distribution and 
some species have become rare and endangered in recent decades 
[12]. Mycorrhizal association is known to be important to orchids 
because they depend on the presence of suitable fungal partners 



for seed germination and seedling development [13,14]. There- 
fore, a complete understanding of the mycorrhizal fungi of the 
many threatened orchid species is required for conservation action 
plans [15]. 

Although molecular taxonomic identification of the endophytic 
fungi of orchid species has now revealed that the diversity of 
orchid associates is much complex [16], the study of the earliest- 
diverging orchid lineages and distribution of mycorrhizal fungal 
associates across orchid phylogeny supported that the ancestral 
state is an association to the Rhizoctonia-like fungi lineages [17], 
and orchid mycorrhizas are predominantly represented by 
associations between photosynthetic plants and Rhizoctonia-Uke 
fungi [15]. The Rhizoctonia-like fungi includes members of the 
Ceratobasidiaceae, Sebacinales and TulasneUaceae [15,17]. How- 
ever, the taxonomy, diversity and distribution patterns of 
mycorrhizal fungi may vary considerably among different orchid 
species or different populations within the same species, ranging 
from very narrow specificity to little specificity [18-26]. Therefore, 
knowing the identity of the fungi that form mycorrhizas with 
orchids and the specificity of the relationships are of crucial 
importance for orchid ecology and conservation [27]. 
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Figure 1. Locations of sampling sites of Liparis japonica in Nortfieast Ciiina. 

doi:l 0.1 371/journal.pone.01 05573.g001 



Liparis japonica, which is a terrestrial photosynthetic orchid 
species, has a wide geographic distribution in the tropical and 
subtropical areas such as China, Korea, Japan and Russia [28]. 
Historical record indicated that this species had a relatively wide 
distribution in Northeast China. Unfortunately, the number of L. 
japonica populations has declined sharply during the last several 
decades and many populations recorded previously has disap- 
peared [29] . In a previous study, we assessed the genetic diversity 
and population genetic structure of L. japonica in Northeast 
China, and found potential restricted gene flow existed among 
populations due to habitat destruction and fragmentation [29] . To 
conserve surviving wild populations and reintroduce plants into 
declining populations requires a more complete understanding on 
the mycorrhizal fungal partners oi L. japonica species. 

In this study, we investigated in detail the taxonomy and 
diversity of culturable mycorrhizal fungi associated with L. 
japonica in Northeast China using morphological and ITS 
phylogenetic analysis. Seven different populations of L. japonica 



species were examined in order to evaluate how fungi diversity 
varied within and between habitats. Using a culture-dependent 

method, the fungal isolates obtained will be important sources for 
deep study and further utilizing in orchid conservation. 

Materials and Methods 

Study species and sampling 

185 individuals of L. japonica belonging to seven natural 
populations were collected in Northeast China at the flowering 
stage in June 2010 and July 201 1, and the distances between these 
populations varied from 25 to 460 km. Details of locations and 
number of sampled individuals of each population were listed in 
Table 1 and Fig. 1. Plants were excavated with a clod of their 
surrounding soil and transferred to the laboratory. The field 
studies did not involve endangered or protected species, and no 
specific permissions were required for these locations. 
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Table 1. Location of the seven populations of Liparis japonica and the corresponding fungal isolates. 



Sample No. of 



Population Code 


Location 


Longitude (E) 


Latitude (N) 


Altitude (m) 


size 


fungal isolates 


DL 


Dongling park, Shenyang, 
Liaoning 


123"34'39" 


41=49'53" 


81 


51 


25 


GG 


Gaoguan, Benxi, 
Liaoning 


124"02'32" 


41=20' 19" 


274 


30 


30 


GMS 


GuanmGnshan forest park, Benxi, 
Liaoning 


1 ZH- D\J 


41 °07'41 " 


438 


18 


27 


CBS 


Changbaishan nature reserve, 
Jilin 


127'47'55" 


42=30'35" 


741 


10 


18 


LTDZ 


Laotudingzi nature reserve, 
Benxi, Liaoning 


124 53'56" 


4r 18'09" 


817 


24 


21 


DQG 


Daqinggou nature reserve, Fuxin, 
Liaoning 


122 12'03" 


42=43'55" 


194 


26 


37 


QS 


Qianshan nature reserve, 
Anshan, Liaoning 


123"07'48" 


41 01 '21" 


163 


26 


43 


Total 










185 


201 



doi:l 0.1 371 /journal.pone.Ol 05573.t001 



Fungal isolation and nnorphological identification temperate, terrestrial orchids wliich have tliick roots and often 

Potential mycorrhizal fungi were isolated from the orchid plants produce abundant pelotons, has few active pelotons suitable for 
and identified from pure culture. L. japonica, unlike many other isolation. Thus, endophytic fungi were isolated from single hyphal 




Figure 2. The nine morphotypes of endophytes isolated from Liparis japonica at seven populations, (morphotype I: Rhizoctonia sp., 
morphotype II: Phomopsis sp., morphotype III: Verticillium sp., morphotype IV: Fusarium sp., morphotype V: Chaetomium sp., morphotype VI: 
Gliodadium sp., morphotype VII: Cylindrocarpon sp., morphotype VIII: Phialophora sp., morphotype IX: Paecilomyces sp.). 
doi:1 0.1 371/journal.pone.01 05573.g002 
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tips emerging from sterilized root portions as in the isolation 
method described by [30]. Three roots per plant were carefully 
cleaned from the soU under running water, surface-sterilized rn 
0.1% (v/w) Mercuric chloride and 75% [v/v) ethanol for 3 min 
and 5 s respectively, and subsequently washed three times in sterile 
distilled water. Root sections of 3-5 mm thickness were obtained 
by cutting and placed in a petri dish with potato dextrose agar 
(PDA). Petri dishes were incubated at 25°C in the dark and 
observed for fungi growing every 2 days for at least 3 weeks. The 
growing colonies were separated onto fresh media for purity and 
this process was repeated three times. Classification of the 
endophytic fungi was based on their growth rate and morpholog- 
ical characteristics, including colonial morphology, production of 
conidiogenous cells, conidial size and dimension on PDA medium 
[31], and similar isolates were grouped into one morphotype. The 
Rhizoctonia-Vike fungal endophytes were recognized by the 
following characteristics: hyphas hyaline with constricted branch 
points, 2.5-7 |j,m diam; submerged growth in PDA; ellipsoid, 
globose or irregular monUioid cells 10-20x10-20 |lm; colony 
creamy white to pale tan or orange, rubbery or leathery in 
appearance and texture [32-34]. 

DNA extraction, amplification and sequencing 

Genomic DNA was extracted from each Rhizoctonia-like fungi 
isolate as described by [35]. ITS sequences of rDNA were 
amplified using the following primer pairs: ITS1-OF/ITS4-OF 
[36]. PCRs were performed in 50-|J,L volumes containing 
approximately 10 |a,g of fungal genomic DNA, 1 xPCR buffer 
(Sangon Biotech Co. Ltd), 1.5 mM MgClz, 200 nM each of dNTP 
(dATP, dCTP, dGTP, and dTTP), 100 nM of each primer and 
2.5 unit ofTaq DNA polymerase (Sangon Biotech Co. Ltd). PCRs 
were carried out with a BIO-RAD T 1 OO"" Thermalcycler with an 
initial incubation at 95°C for 5 min, followed by 35 cycles with 
denaturation at 95°C for 30 s, annealing at 60°C for 1 min and 
extension at 72°C for 1 min, and a final extension reaction at 
72°C for 10 min. The PCR products were purified using the 
TaKaRa Agarose Gel DNA Purification Kit (TaKaRa, Dalian, 
China) and sequenced using the BigDye terminator v3.1 on an 
ABI 3730x1 DNA Analyzer (Applied Biosystems, Forster City, 
California, USA). New ITS sequences from this study have been 
deposited in GenBank (the accession numbers is from KF537635 
to KF537658). 

Data analysis 

ITS sequences were analyzed with BLAST [37] against the 
NCBI sequence database to find the closest sequence matches in 
the GenBank database (http://www.ncbi.nlm.nih.gov). Genus and 
species of the database match were accepted whenever identity 
between our sequence and that of the database was greater than 
97% [38]. Sequences used to construct phylogenetic trees were 
from organisms considered to be closely related to the mycorrhizal 
fungus of L. japonica based on BLAST searches, and additional 
orchid mycorrhizal fungi sequences (mainly based on [19]) were 
also chosen for phylogenetic analysis. Selected sequences were 
aligned using the program Clustal_X 2.1 [39]. Alignments were 
checked for ambiguities and adjusted manually when necessary 
and the alignment gaps were treated as missing information. 

Phylogenetic analyses were estimated using maximum parsimo- 
ny (MP), Neighbor-joining (NJ) and Bayesian inference (BI) 
methods in PAUP* version 4.0bl0 [40], MEGA version 4 [41] 
and MrBayes version 3.0b4 [42], respectively. For MP analysis, 
the bootstrap support values were estimated by 1000 replicates 
with 10 random sequence additions and tree-bisection-reconnec- 
tion (TBR) branch swapping. For NJ analysis, genetic distances 
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Table 3. Morphological characters of the 9 morphotypes of fungal isolates from Llparls japonica. 





Morphotype 


Colony color 


Colony 
texture 


Growth rate 
(mm day^) 


Aerial mycelium 


Conidial 
shape 


Conidia 
Size [\im) 


1 


off-white 


leathery, compact 


11.9±1.5 


absent 






II 


pink to tan 


leathery, 

radial 

draped 


4.9±0.5 


absent 


elliptic cylindrical 


6-8x3-4 


III 


white or buff 


cottony, 
slightly loose 


5.7±0.4 


scarce 


fusiform 


5x1-1.5 


IV 


white 


lanose, 
loose 


11.8±0.5 


abundant 


ovate, 
cylindrical 


cylindrical: 

15-20x5-10 
ovale: 4x3 


V 


white at 
first, 

then olive-green 


pulvinate, compact 


4.8±0.5 


absent 


llmoniform 


8-9 x4-5 


VI 


white 


cottony, 

slightly 

loose 


5.0±0.2 


abundant 


elliptic 


6-8 x3-5 


VII 


pale tan 


cottony, compact 


5.2±0.2 


abundant 


clavate, 
elliptic 


40x5-10 


VIII 


Pale taupe 


felted, 
compact 


3.4±0.4 


scarce 


ovate 


3-5x2.5 


IX 


lilac 


felted, 
compact 


3.3±1.7 


abundant 


elliptic 


4-5x2 



doi:10.1371/journal.pone.0105573.ta03 



were calculated using the Kimura two-parameter model, and the 
bootstrap support values were estimated by 1000 replicates. Prior 
to the Bayesian analysis, the Akaike information criterion (AIC) in 
Modeltest version 3.7 [43] was used to select the best-fit model of 
molecular evolution for each dataset. For BI analysis, four chains 
of the Markov Chain Monte Carlo (MCMC) chains were run, 
sampling one tree every 1000 generations for 3,000,000 genera- 
tions, starting with a random tree. Bayesian posterior probabilities 
(PP) were calculated from the majority rule consensus of the tree. 
The sequence of Septobasidium carestianum (GenBank accession 
number DQ241448) was selected as outgroup [21]. 

Results 

Fungal distribution and morphological diversity 

After isolation and purification, a total of 201 isolates of 
endophytic fungi from 185 plants were obtained from seven L. 
japonica populations. According to their morphological characters 
and growth rate on PDA medium, fungal isolates were classified 
into nine morphotypes with four to seven morphotypes in each 
population (Table 2). The isolates in morphotypes I and II were 
observed in all populations studied and represented the majority of 
the fungal community (61.7%), while isolates from other 
morphotypes were not common and occurred only at one or a 
few populations (Table 2). Morphological characters and detailed 
descriptions of the nine morphotypes were given in Fig. 2 and 
Table 3, and according to the colony and micro-morphological 
characters, we identified morphotype I as Rhizoctonia sp., 
morphotype II as Phomopsis sp., and morphotype III to IX as 
Verticillium sp., Fusarium sp., Chaetomium sp., Gliocladium sp., 
Cylindrocarpon sp., Phialophom sp., and Paecilomyces sp., 
respectively. Only the Rhizoctonia-like isolates in morphotype I 
were given codes and subjected to further collection of mycelium 
for DNA extraction and phylogenetic analysis. 

Of the 185 individuals of L. japonica samples, each of 180 
samples gave rise to one Rhizoctonia-]ike isolate, whereas the other 



five samples gave rise to two Rhizoctonia-like isolates, and no 
sample gave rise to three or more Rhizoctonia-Vikt isolates. 

Phylogenetic analysis 

About 640-bp portions of the ITS were sequenced for the 24 
representative isolates of morphotype I, with 3 to 4 isolates within 
each population. All sequences identified by BLAST in NCBI 
were from Tulasnellaceae strains. Thus, the representative 
Tulasnella spp. that have been identified were selected for 
phylogenetic analysis with the sequences obtained in this study, 
and a total of 54 different sequences were used. 

The overall topology of the trees produced with Maximum 
parsimony (MP), Bayesian inference (BI) and Neighbour-joining 
(NJ) methods were similar based on our ITS data (Fig. 3). The MP 
strict consensus tree (Fig. 3) showed that fungal associates of L. 
japonica all belonged to the genus Tulasnella, and all isolates 
belonged to a single major clade. Some mycorrhizal isolates from 
other orchid species which have been identified as T. calospora, 
and a mycorrhizal isolate {Epulorhiza sp., AM040890; [44]) from 
L. loeselii species in Hungary also fell in this major clade. 

Discussion 

Compared with the direct DNA amplification technique, 
identification of orchid mycorrhizal fungi with the pure culture 
technique has been criticized because of the bias toward fast- 
growing and easily cultured fungi [15]. However, the culture- 
dependent technique is still used in a vast majority of orchid 
mycorrhizal fungi studies, mainly because fungi partners of these 
plants were usually easy to isolate and no significant differences 
were found between the pure culture technique and direct DNA 
amplification technique [20,23,45-47]. Unlike many other tem- 
perate, terrestrial orchids which have thick roots and often 
produce abundant pelotons, L. japonica has few active pelotons 
suitable for isolation and direct PCR identification, so we used 
standard culture method for fungi identification in this study, and 
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72/0.95/- 
72/0.96/- 
-/-/- 



-/-/- 



90/-/- 



-/0.5/- 



100/0.93/100 



100/1/100 



81/0.52/73 



-/0.52/58 



55/0.57/62 



-/0.76/- 



54/-/ 



83/1/69 1 
99/0.99/81 r 



77/0.54/ 



-/-/ 



100/1/99 



65/1/64 



98/1/100 



97/0.72/77 r 
100/1/991 L 



-tv 



97/ 0.72/77 r 
100/1/99 L 



93/0.98/- 



IOCn/100 



100/1/100 



LTDZ 8 
DQG 16 
DL 52 
LTDZ 6 
LTDZ 14 
LTDZ 2 
GMS 87 
QMS 99 
QMS 67 
GMS 65 
CBS 79 
CBS 104 
GG 75 
CBS 86 
CBS 114 
QS 77 
QS 80 
GG 69 
DQG 101 
DQG 29 
DQG 19 
GG 106 
GG71 
QS 91 

AM040890 Epulirhiza sp.* 
GU166403 Tulasnella calospora 
Fjei3262 Tulasnella calospora 
GUI 66407 Tulasnella calospora 
EF393622 Tulasnella calospora 
AY373291 Tulasnella calospora 
EU218888 Tulasnella calospora 
AY643804 Tulasnella calospora 
AY373290 Tulasnella bifrons 
DQ388044 Tulasnella calospora 
AY373288 Tulasnella sp.* 
AY373287 Tulasnella sp.* 
DQ388041 Tulasnella calospora 
DQ388045 Tulasnella calospora 
EF393631 TulasrmOa calospora 
AY373298 Tulasnella calospora 
AY373297 Tulasnella danica 
EU2 18891 Epulorhiza anaticula 
AY373292 Tulasnella eichleriana 
AY373293 Tulasnella violea 
AY373296 Tulasnella tomaculum 
EU218889 Tulasnella irregularis 
AY373303 Tulasnella violea 
DQ457643 Tulasnella violea 
AY373294 Tulasnella albida 
DQ457642 Tulasnella pruinosa 

' AY373295 Tulasnella pruinosa 
DQ388048 Tulasnella asymmetrica 

- DQ^8047 Tulasnella asymmetrica 
DQ241448 Septobasldium carestianum 



Figure 3. Strict consensus tree of Maximum parsimony (iVIP) based on ITS sequences of representative isolates of Rhizoctonia-Wke 
fungi isolated from Liparis japonica and identified Tulasnella spp. and Epulorhiza spp. Numbers at the nodes are MP bootstrap percentages, 
Bayesian posterior probabilities and Neighbour-joining support percentages (£50%), respectively. "-" indicates node not supported in Bayesian 
inference and Neighbour-joining analyses. An asterisk indicates fungi isolated from Liparis lilifolia and L. loeselii (other Liparis species). The 
Septobasldium carestianum sequence was designated as outgroup for rooting the trees. Sequences obtained from GenBank were shown with 
accession numbers. 
doi:1 0.1 371/Journal.pone.01 05573.g003 
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the results were credible by the striking observation of the same 
fungal clade isolated across such a wide geographical distribution 
of L. japonica. Moreover, although the molecular method is 
powerful for fungi identification, isolation and cultivation of the 
mycorrhizal fungi isolates are still needed for deep study and 
further utilizing in orchids conservation. 

In the present study, various endophytic fungi were isolated 
from roots of L. japonica and identified by morphological and the 
rDNA ITS sequences, and most of the identified endophytes of L. 
japonica are Basidiomycetes and Ascomycetes. In general, 
endophytic fungi in orchid plants are ver\' abundant, but it seems 
likely that only RhizoctoniaASkc fungi species may form mycor- 
rhizal associations with host [48] . In this study, 86 endophyte fungi 
isolates belonging to the Tulasnella genus were isolated (42.8%) 
based on morphological characters and phylogenetic relationships 
recc)nstruct(;d from ITS sequence, and they existed in every L. 
japonica population. Although some endophytes have been shown 
to be able to increase the tolerance of some plant species [49], little 
is known about the role of endophytes for orchid performance, 
and the non-RhizoctoniaAike. endophytic fungi also deserve 
additional study. 

The specificity of orchid mycorrhizal associations has important 
implications for orchid biolog}', conservation and restoration of 
orchid populations [50]. Specific fungal partners may lead to 
enhanced seed germination rates, and the cfri( i(;ncy of nutrient 
exchange between partners may be heightened with specific plant- 
fungus combinations [45]. Our results showed that all mycorrhizal 
fungi isolates obtained belonged to a single clade of Tulasnella, 
and the molecular analysis showed that they are all representatives 
of the same species, T. calospora, which indicates that the 
relationship between L. japonica and its mycorrhizal fungus is very 
specific regardless of habitats, at least in Northeast China. There 
are more and more reports showed that the degree of specificity 
can vary among species independent on orchid types [19,22]. 
There is also evidence that a preference exists among some orchids 
for specific fungal partners, and some photosynthetic orchids, even 
when sampled over a wide range, have a single dominant 
mycorrhizal fungus, for example, Spiranthes sinensis [33], Good- 
yera pubescens [51], and Pterostylis nutans [52]. Our results also 
supported that high specificity in widespread and photosynthetic 
orchids may be more common than previously thought. More- 
over, the phenomenon that multiple fungal taxa existed in a single 
root system was detected in some photosynthetic orchids which 
indicated that multiple fungi in a single root system may occur 



more often in photosynthetic than in mycoheterotrophic orchids 
[26]. While our results showed that this characteristic is not 
universal in orchid mycorrhizal systems. 

From another perspective, the mycorrhizal specificity we 
observed might not be accurate for the reason that our results 
were based on a single collection of root tissue in adult L. japonica 
plants. However, the diversity and identify of orchid mycorrhizal 
fungi may vary according to different life stages of plant. 
Specificity during germination and early development has been 
detected in orchids, whereas the adult stage of plant growth in the 
same species can have less specific fungal partners [13,51,53]. 
Accordingly, fungi required for germination and recruitment of L. 
japonica might be different from those associated with adult plants, 
and this will be the next logical steps to investigate the pattern of 
mycorrhizal association. 

Although L. japonica had a strong preference for T. calospora 
species across the sampled range, the fungus did not have an 
equivalent specificity. T. calospora has been shown to be a very 
common mycorrhizal fungus of many terrestrial orchid species 
[15], and the fact that T. calospora stains were present in different 
populations of L. japonica in Northeast China in our study and 
many other orchids in other regions of the world indicates its wide 
distribution in different ecological environments. Because of this, 
we considered that the decline of L. japonica in Northeast China 
was not resulted from fungal specificity, although there exists a 
point that a dependence (m narrowly specific interactions with 
fungi and pollinator may predispose many orchids to become rare 
[15,45,54]. 
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